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Abstract 

The benchmark problems in Category 1 (Internal 
Propagation) of the third Computational Aeroacous- 
tics (CAA) Workshop sponsored by NASA Glenn 
Research Center are solved using the space-time 
conservation element and solution element (CE/SE) 
method. The first problem addresses the propaga- 
tion of sound waves through a nearly choked tran- 
sonic nozzle. The second one concerns shock-sound 
interaction in a supersonic nozzle. A quasi 1-D 
CE/SE Euler solver for a nonuniform mesh is de- 
veloped and employed to solve both problems. Nu- 
merical solutions are compared with the analytical 
solution for both problems. It is demonstrated that 
the CE/SE method is capable of solving aeroacous- 
tic problems with/without shock waves in a simple 
way. Furthermore, the simple non-reflecting bound- 
ary condition used in the CE/SE method which is 
not based on the characteristic theory works very 
well. 

1. Introduction 

The method of space-time conservation element 
and solution element (abbreviated as the CE/SE 
method) is an innovative numerical method for solv- 
ing conservation laws. It is different in both con- 
cept and methodology from the well-established tra- 
ditional methods such as the finite difference, finite 
volume, finite element and spectral methods. It is 
designed from a physicist’s perspective to overcome 
several key limitations of the traditional numerical 
methods. 

Simplicity, generality and accuracy are weighted 
in the development of this method while the funda- 
mental requirements are satisfied by the scheme. Its 
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salient properties are summarized briefly as follows. 
First, the concepts of conservation element and solu- 
tion element are introduced to enforce both local and 
global flux conservations in space and time instead 
of in space only. Second, all the dependent variables 
and their spatial derivatives are considered as indi- 
vidual unknowns to be solved for simultaneously at 
each grid point. Third, no approximation techniques 
other than Taylor’s series expansion, no monotonic- 
ity constraints, and no characteristic-based tech- 
niques are used in the design of the scheme. A de- 
tailed description of this method and the accompa- 
nied analysis can be found in [1-3]. 

A variety of numerical tests have been performed 
previously to illustrate the accuracy of this method. 
For the CE/SE Euler solver, highly accurate numer- 
ical solutions have been obtained for various flow 
problems involving discontinuities, such as shock 
waves, contact surfaces and even their interac- 
tions^]. Moreover, applications of the same Euler 
solver to computational aeroacoustics (CAA) prob- 
lems reveal that the accuracy of the results is com- 
parable to that of a 4th-order compact difference 
scheme even though the current solver is only 2nd- 
order accurate, and the non reflecting boundary con- 
dition can be implemented in a simple way without 
involving characteristic variables. It can be applied 
to subsonic, transonic, and supersonic flows in the 
same form without using the characteristic-based 
techniques. Results show' that the present solver can 
handle both continuous and discontinuous flows very 
well[5 -1 1] . 

In this paper, a quasi 1-D CE/SE Euler solver 
for a nonuniform mesh is developed and used to 
solve the benchmark problem in Category 1 (Inter- 
nal Propagation) of the third Computational Aeroa- 
coustics (CAA) Workshop. Among the two proposed 
problems, one addresses the propagation of sound 
weaves through a transonic nozzle and another con- 
cerns shock-sound interaction in a supersonic noz- 
zle. For problem 1, both uniform meshes with con- 
stant numerical dissipation models and nonuniform 
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meshes with variable numerical dissipation models 
are used. Numerical results are presented and com- 
pared with the analytical solution. A good agree- 
ment is achieved by using a stretched 401 point 
mesh. For problem 2, 101 and 201 point uniform 
meshes are used and excellent agreements with the 
analytical solution are obtained. Further, the nu- 
merical solutions obtained by solving the linearized 
equations are compared with those obtained by solv- 
ing the nonlinear equations for both problems. 


u, respectively, being the area and the outward unit 
normal of a surface element on 5(1'); and (iv) dv is 
the volume of the space-time region V. Note that (i) 
hm ’ds is the space-time flux of h m leaving the region 
V through the surface element ds\ (ii) f f v g m dv is 
the volume integration of the source term g m over 
the space-time region V\ and (iii) all mathematical 
operations can be carried out as though E 2 were an 
ordinary two-dimensional Euclidean space. 

As a preliminary, let 


2. Governing Equations and Numerical Schemes 

In the following, the nonlinear quasi ID Euler 
equations are described in conservative form first, 
which is followed by numerical schemes, numerical 
results and conclusions. 


fm,i = dfm/dui, m,l= 1,2,3. (10) 

The Jacobian matrix F, which is formed by 
fm,h J = 1, 2, 3, is referred to in [1], 

Because f, n . m = 1,2,3 are homogeneous func- 
tions of degree one in 10,7/2, and 7/3, we have 


2.1 Governing Equations 


The quasi 1-D Euler equations in conservative 
form are expressed as 


Olim & fm 

IT + ~W ~ 9m 


m = 1,2,3 


(i) 


in which 


Ui-p, u -2 = pv, u 3 = p/(j - 1) + pir/2 (2) 


arid 


/l = U -2 (3) 

h = (7 ~ 1)«3 + (3 - 7)(w 2 ) 2 /(2wi) (4) 

h = 7«2«3/«i - (l/2)(7 - 1 )(u 2 ) 3 /(ui) 2 (5) 

dA 1 

91 = -“•te 3 (0) 

= P) 

, dA 1 

93 ~ dx A ^ 


where p, t\p, and 7 are the mass density, velocity, 
static pressure, and constant specific heat ratio, re- 
spectively, and A is the area of a nozzle. 

Let xi — x and X 2 = t be the coordinates of a 
two-dimensional Euclidean space, F2. The integral 
form of Eq. (1) in space- time E 2 is 


/ 

JS(V) 



777 = 1,2,3 


(9) 


where (i) S( V) is the boundary of an arbitrary space- 
time region V in F 2 ; (ii) h m = (/ m , n m ), m = 1, 2, 3, 
are the space-time mass, momentum, and energy 
flux vectors, respectively; (iii) ds = da n with da and 


3 

fm — ^ ] fm,l u l- 
7=1 


( 11 ) 


2.2 The ID Space-Time Nonuniform Mesh 

Consider a nonuniform mesh shown in Fig. 1, 
in which the solid circles are the grid points and 
crosses are the solution points. Assume the spa- 
tial coordinate of the grid points be xj with j — 
0, 1/2, 1, 3/2, 2, 5/2. Generally, the grid point j is 
not the centroid of the line formed by grid points 
j -h 1/2 and j — 1/2. Let x'- denote the spatial coor- 
dinate of the solution points. Define 

x j = ( x j- 1/2 + * j+ i/a)/2 (12) 

for interior mesh points, while x' = xj for boundary 
points. 

The definitions of the conservation element (CE) 
and solution element (SE) at a grid point j are the 
same as those for ID uniform mesh in [1]. The 
flow marching variables are computed at the solution 
points, while the grid points are used to construct 
CEs and SEs. For later use, we define 


. _ def 

Axj = xj - Xj_ l/2 

(13) 

Ax+ d = Xj+I/ 2 - Xj 

(14) 

A def 

LlXj — Xj+1/ 2 F/ — 1/2* 

(15) 


2.3 The Quasi ID Euler a Scheme 

For ^ any (x,t) € SE(j,n), u m (x, f), / m (x,t), 

and hm(x, t) are approximated by 
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fm(x, t‘,j,n), and h* m (x,t\j,n), respectively. They 
will be defined shortly. Let 


UmOM; j,n) = + + 

(16) 

where (u m )*?, (u m;r )", and (u m f)J are constant 
within SE(j, n). Obviously, they can be consid- 
ered as the numerical analogues of the values of 
u mj du m /dx, and du m /dt at (xj,t n ), respectively. 

Let {f m )j and denote the values of f m and 

f m j , respectively, when u m ,m = 1, 2, 3, respectively, 
assume the values of (u m )",m = 1,2,3. For any 
?7i = 1, 2, 3, let 


and 


Because 


and 


W)" = £(/m,i)?(«,x)7 
1=1 

(17) 

3 

f,n<)] = 'EifmtfM] 

(18) 

1=1 


df m ^ j. 

a at 

(19) 

£=1 


dfm ^ 

/-I 

(20) 


(frnxYj and ( fmt)j c an be considered as the numer- 
ical analogues of the values of df m /dx , and df m /dt 
at(xj,t n ), respectively. As a result, we define 


/* (x.ijj.n) d = (/ m )?+(/ m x)J(x-x')+(/ mt )J(#-< n ) 

( 21 ) 

Because ft m = (/ m ,u m ), we also assume that 

h* m (x 9 t;j,n) = (f^{x y t;j,n),u; n (x y t;j,n)) (22) 

Note that, by their definitions: (i) (f m )j » 
are functions of (u m )”, m = 1,2,3; (ii) (fmx)j are 
functions of (u m )” and (u mx )™, m = 1,2,3; and 
(iii) (/mOj are functions of (ti m )” and (u m ,)J, m = 
1,2,3. 

Moreover, we assume that for any (:r, i) € SE(j, n), 
and for m = 1, 2, 3, 

du* m (x,t;j,n) , df^{x,t;j,n) , /0<n 

dt + & ~ ( ^ 

Note that Eq. (23) is the numerical analogue of 
Eq. (1). With the aid of Eqs. (16) and (21), Eq. 
(23) implies that, for any m = 1,2,3, 

(' u mt)j “ ~(/mx)j + (#?n)j 

= -E[/^««»]; + (flm)y ( 24 ) 
?=1 


Since (^ m )” is also a function of (u m )", m - 1,2,3, 
thus (u m t)j is a function of (w m )" and ( u mx )j . From 
this result and the facts stated following Eq. (22), 
one concludes that, for each m — 1,2,3, there are 
only two independent discrete variables, i.e., (u m )", 
(u tni )J, associated with the mesh point (j,n). 

Further, the flux h* n is conserved over CE ± (j,n), 
i.e., 


l 


S(CE- (j.n)) 


K> 


■ ds = (g m )"_ Ijl + (g m )j 


i A*A*J 


(25) 


and 


AtAxj 


<£ h m ■ ds = \(g m ) j+1 f 2 + (g„ 

JS(CE+U, n)) L 

(2G) 

where S(CE±(j,n)) denotes the boundary of 
CE±(j,n). Note that Eqs. (25) and (26) are the 
numerical analogues of Eq. (9) . 

Let Uj,(u x )j and gj , respectively, be the 3x1 
column matrices formed by (u m )”, (u mx )j and 
(^m)j , m = 1,2,3. Let Fj 7 denote the 3x3 Jacobian 


matrices formed by (/ m ,j)j, m,/ = 1,2,3. Define 


(f + )j = f ^f; 


(27) 


(^)? d = ( 28 ) 

With the above preparations, a new set of 3 x 3 
matrices will be defined at each mesh point (j, n). 
To simplify notation, we will strip from all 3 x 3 
matrices involved, their indices j and n. Define 


Sfi 

def 

Axj I + 

2F + 



(29) 

Sfa 

def 

Axj(^ 

-1/2 + 
2 

fl 





2 (F + ) 2 




(30) 


def 

Ax"/ + 

2F + 



(31) 

Si2 


l> 

i 

I* 

-1/2 + 
2 


~ x j— 1 / 2 )^ _ 




2(F+) 2 




(32) 

S 2 + ! 

def 

Ax+I- 

2F + 



(33) 


def 

Ax+C J+1/2 4 

J v 2 

_ 

-x))I- 




2(F+) 2 




(34) 

^21 

def 

A art/- 

2F+ 



(35) 

^22 

def 


fl . 

~ x 'j+ 1 / 2) 1 + 




2(F+) 2 




(36) 
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where I is the 3x3 identity matrix. 

Note that, with the aid of Eqs. (12)— (15), it can 
be shown that 


the 2D CE/SE Euler a scheme is formed by Eq. (41) 
and 


(**)? - W)] 


(45) 


s n + £« = (37) 

E+ + E+ = 0 (38) 

Next we will show how to compute Eqs. (25) and 
(26). Because u* n (x,t;j,n) is linear in x and £, the 
total flux of h* m leaving CE±(j, n) through any one 
of the line segments that form its boundary is equal 
to the scalar product of the vector h* m evaluated at 
the centroid of the line, and the normal vector of the 
line. The results of this evaluation can be combined 
into the matrix form 


[E+j (u - g*) + Et 2 «x]" = 

[Efi (<? + t) + (39) 

and 

[E+(u-^) + E+^]; - 

[s 2 - I («+r)+s 2 - 2 ^];; 1 1 / / 2 2 . (40) 

By adding the two matrix equations and using 
Eqs. (37) and (38), one concludes that 

= £n (* + 1) + 

(41) 

[E5i«r+r)+2itf«];; a l /2 2 } 

i.e., Hj can be obtained using an iterative method, 
such as Newton iterative method, in terms of the 
marching variables at the (n — l/2)th time level. 
Furthermore, it can be seen that the matrices 
(E+)”,r, s = 1,2 and are functions of Uj 

only, they become known once Uj is known. Thus 
{u x )j can be computed using either Eq. (39) or (40). 
Consider Eq. (39) for an example. Define 


2.4 The Quasi ID Euler a-e-a-fi Scheme 

The quasi ID CE/SE Euler a-e-a-P scheme for a 
nonuniform mesh can be derived in the same way 
as that for ID CE/SE Euler a-e-a-P scheme for a 
uniform mesh. 

Let 

(0; ±1/ 2 = M n j±l/l + (Ao/2(u ra< );-$ (46) 

and define 

=' [(«m)7 - (0"-,/ 2 ] /AX- 

«nx+)j‘ = [«)] + l/2 ~ («m)j B ] /Ax+ 

and 

«x)? = f [«,,)] +l/2 - (u' m )?-i/a] /A Xj . 

Further, we can define 

. u. = l«mz-l°«m»+ + l M mx+i a »mx- 
mT |ti^_|“ + K, + |« 

for \u% x _\ a + \u c mx+ \ a > 0, otherwise u” T = 0. 

The quasi ID CE/SE Euler a-e-a-P scheme is 
formed by Eq. (41) and 

(**)? = (^)"+2e((^)"-(^)p+i3((t^)^-(^)7) 

(51) 

where e ) a and (i are parameters for controlling nu- 
merical dissipations. They can be either constants 
or variables. With e, a and p being held constant, 
numerical dissipation associated with the present 
scheme at a mesh point (j, n) tends to increase as 
the local Courant number i/* decreases. To com- 
pensate for this effect, e and ft in Eq. (51) may be 
replaced by e(v'j') and /3(i/”) which are expressed as 


(47) 

(48) 

(49) 

(50) 


(&)7 = [Sfi(tf + t) + - 1 )]; ■ 

(42) 

Then Eq. (39) can be rewritten as 

[Sf 2 i?j]" = (Si) j. (43) 

With assumption that [(Ef 2 )”] 1 exists, we can 
have 

(tf x )7 = [(E+)7] -1 (^)y. (44) 

Note that ('«*)" can be computed in the similar way 
using Eq. (40). Let («“)" denote the solution. Then 


e({/") = 0.5i/"exp(l — i /") and 

0W) = yft, 0 < v? < 1. (52) 

Note that e(v™) and /3(i/7), are monotonically in- 
creasing function of i/j with c(0) = 0 and /3(0) = 0. 
A ID CE/SE Euler solver with and P{Vj) in 
Eq. (52) was used in [1] to obtain accurate shock 
tube solutions with the maximal Courant number 
in each numerical simulation ranging from 0.88 to 
0 . 022 . 

Similarly, the quasi ID CE/SE Euler solver is 
second-order accurate in both space and time. The 
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stability condition is the Courant number less than 
unity. 

3. Numerical Results 

For the two problems solved in the following, the 
CE/SE quasi ID Euler solver described previously 
for the nonlinear equations is used to compute the 
steady-state solution of the nozzle flowfield. For the 
simulation of acoustic wave propagation, both the 
linear and nonlinear quasi ID Euler solvers are used 
under the same computational conditions. The lin- 
ear solver is constructed for the linearized quasi ID 
Euler equations in conservative form. It has the 
identical formulation as the nonlinear solver except 
the time-marching variables are different. The de- 
tails of the linear solver will be described in a sepa- 
rate paper. Both the linear and nonlinear numer- 
ical solutions are presented. Numerical solutions 
of the acoustic wave plotted in all figures are non- 
dimensionalized by its amplitude £. In the computa- 
tions using a uniform mesh, the constant numerical 
dissipation model with e = 0.2, a = 0, and ft = 1 is 
used for both problems. 

3.1. Propagation of Sound Waves through 
a Transonic Nozzle 

The first problem addresses propagation of sound 
waves in a transonic nozzle, in which the local Mach 
number of flow near the throat maj r be close to sonic. 
The computation of sound propagating through such 
regions presents a challenging problem. The area of 
the nozzle is 

A(x) - 

f 0.536572- 0.198086 exp (-(In 2)(^) 2 ) , x>0 
\ 1.0 — 0.661514 exp (-(In 2)(^g) 2 ) , x < 0. 

(53) 

In the uniform region downstream of the throat, the 
diameter of the nozzle is D and the main flow is 
uniform with Mach number, Moo? of 0.4, speed of 
sound, a OCi and static density, Poq. Flow variables 
are non-dimensionalized by using as the velocity 
scale, D as the length scale, D/a oo as the time scale, 
Poc as the density scale, and p 0 c a4 as the pressure 
scale. Thus the mean flow in the uniform region 
downstream of the throat is 

p — 1, v = 0.4, p= 1/1.4. (54) 

The acoustic wave, with angular frequency uj — 0.6tt, 
is described as 

p'=p , = £C0sL(j-4 i j-+t) , 


X 

v 1 — — s cos u>(~ — — — M) (55) 

1 — nice 

where e = 10“ 4 . It is generated downstream arid 
propagates upstream through the narrow passage of 
the nozzle throat. The computational domain is 
-10 < x < 10. It is recommended that computa- 
tions use no more than 400 mesh intervals. 

As the first step, the steady-state solution of the 
nozzle flow is computed. The initial conditions are 
specified using flow properties in the uniform down- 
stream region. For imposing the boundary condi- 
tions, the back pressure (pback — 1/1.4) is specified 
at the outlet, and the total pressure and total density 
are specified at the inlet. The other needed informa- 
tion at both the inlet and outlet are obtained using 
extrapolation from their neighboring mesh points. 
The steady-state solution obtained using a 401 point 
uniform mesh with CFL = 0.847 is shown in Fig. 2 
in comparison with the exact solution represented 
by solid lines. It can be seen that flow properties 
are uniform in most region of the nozzle, but change 
dramatically near the nozzle throat. 

After the steady-state nozzle flowfield is available, 
the acoustic wave propagation can be simulated us- 
ing the same nonlinear Euler solver. The initial con- 
ditions are specified using the steady-state solution 
of nozzle flowfield. The boundary condition is dif- 
ferent from that used for computing the steady-state 
solution. At the outlet, the time-marching variables 
are specified as 


(Ml)" = P + P' 

(56) 

(u 2 )] = (p + p')(v + v 1 ) 

(57) 

—r + \(p + p')(v + v') 2 . 
7 — 1 2 

(58) 


At the inlet, the non-reflecting boundary condition 
is enforced by setting 

(t/ m )" +1/2 = (u m )” +1/2J (u mi )J +1/2 =0 (59) 

where m = 1, 2, 3. 

Different mesh sizes are tested in the current com- 
putation. First, uniform meshes with Ax = 0.05 and 
0.025 in the computational domain of — 10 < x < 10 
are used. Numerical solutions of the acoustic pres- 
sure p 1 at t = 20T are shown in Figs. 3 and 4, respec- 
tively and compared with the analytical solution. It 
can be seen that the peak values near the throat 
and ■wave pattern upstream of the throat can not be 
captured correctly using a 401 point uniform mesh, 
while a greatly improved result is produced using 
an 801 point uniform mesh. The amplitude of the 
acoustic wave upstream of the throat is lower than 
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that given by the analytical solution. Further, the 
corresponding solutions obtained with Ax = 0.0125 
in the computational domain of -5 < x < 5 are 
shown in Fig. 5. Although the peak values are 
still not the same as the analytical solution, the pro- 
files of the p r both upstream and downstream of the 
throat agree very well with the analytical solution. 

In order to reduce the number of grid points, a 401 
point nonuniform mesh that cluster near the throat 
is used. The ratio of the largest to smallest mesh in- 
terval is around 10 with Ax m in — 0.0047. The vari- 
able €j and /?” described in Eq. (52), and constant 
a — 0 are used in the computation. The obtained 
acoustic pressure at t — 20 T with CFL = 0.9084 
is shown in Fig. 6 which is nearly similar to that 
obtained using an 801 point uniform mesh in the 
domain of — 10 < x < 10. Further improvement can 
be made in the future by reducing the numerical dis- 
sipations or using a more stretched mesh. Its conver- 
gence to the time-period solution is shown in Fig. 7. 
Under the same computational conditions and mesh, 
the corresponding solution obtained using e = 10~ 3 
is plotted in Fig. 8 showing skewness upstream of 
the throat. It can be concluded that the amplitude 
of the acoustic wave should be small enough to avoid 
the nonlinearity when a nonlinear scheme is used to 
obtain the linear solution. In the present computa- 
tion using a nonlinear solver, the steady-state solu- 
tions converge to 10“ 7 , and e = 10 -4 is used. 

For the linear solver, the amplitude of the acous- 
tic wave is set to 1 and the steady-state solution 
obtained by using the nonlinear solver is used in the 
computation. All the three uniform meshes used for 
the nonlinear solver described above are tested with 
the linear solver. The linear solutions are very simi- 
lar to the nonlinear solutions. The solution obtained 
using an 801 point uniform mesh in the domain of 
— 10 < x < 10 is shown in Fig. 9 as an example. The 
nonlinear solution is slightly better than the linear 
solution under the same computational conditions. 


3.2 Shock-Sound Interaction 


In this problem, the same nozzle geometry from 
the previous problem is used, but here there is a su- 
personic shock downstream of the throat. All quan- 
tities are non-dimensionalized using the upstream 
values. The Mach number at the inlet, M 0 CJ is 
0.2006533 and the back pressure, pback » is 0.6071752. 
Following the same steps described in the first prob- 
lem, the steady-state nozzle flowfield is computed 
and used as the initial condition for the sound-shock 
interaction simulation. The acoustic wave is de- 
scribed as 


p — p — v — £ sin 




1 + Mry 


+ 0 


(60) 


where e and w have the same values as that in the 
first problem. At the inlet, the time-marching vari- 
ables are specified using Eqs. (56)-(58), while at the 
outlet, for 77? — 1,2,3, 

(u m )] +1/2 = («m)7- J/2 , («mi)J +l/J = 0 (61) 

are used as the non-reflecting boundary conditions. 

Uniform meshes with 101 and 201 grid points are 
used for this problem. The steady -state solution ob- 
tained using a 101 point uniform mesh is shown in 
Fig. 10 and compared with the analytical solution 
represented by solid lines. Numerical solutions of 
the acoustic wave obtained using 101 and 201 point 
uniform meshes are shown in Fig. 11, in which only 
the analytical solution of p f is plotted for a com- 
parison. The analytical solutions of the density and 
velocity are not available. It can be seen that the so- 
lution obtained on a 101 point uniform mesh is very 
close to the analytical solution and that obtained 
on a 201 point uniform mesh appears to graphically 
match the analytical solution. A very solid conver- 
gence to the time-period solution is demonstrated in 
Fig. 12 for the coarse mesh. Further, the acoustic 
pressure at the outlet in one time period is plotted 
in Fig. 13 along with the analytical solution show- 
ing an excellent agreement. The linear solution is 
almost identical to the nonlinear solution which is 
shown in Fig. 14 for a 101 point uniform mesh. 


4. Conclusions 


The quasi 1-D CE/SE Euler solver for a nonuni- 
form mesh is developed and applied to simulate the 
propagation of sound waves and shock-sound inter- 
action in a nozzle. For problem 1, a satisfactory nu- 
merical result is obtained using a nonuniform mesh 
with 401 grid points. For problem 2, an excellent 
agreement between the numerical results obtained 
on a 101 point uniform mesh and the analytical solu- 
tion is achieved. It is demonstrated that the CE/SE 
method is capable of solving aeroacoustic problems 
with/ without shock waves in a simple way. No spe- 
cial techniques are used for shock waves, and acous- 
tics waves and shock waves can be captured concur- 
rently within the computational domain. 
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t(n) 
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solid line: exact sol. 

delta: 401 point uniform mesh. 



Figure 1: A nonuniform staggered mesh in the spare- 
time domain. 


-4 









Figure 4: Comparisons between the exact solution Figure 5: Comparisons between the exact solution 

and CE/SE solutions obtained using an 801 point and CE/SE solutions obtained using an 801 point 

uniform mesh in -10 < x < 10. uniform mesh in -5 < x < 5. 
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Figure 11: The CE/SE nonlinear solutions of the 
acoustic wave at t = 40T using 101 and 201 point 
uniform meshes. 


Figure 10: The steady-state CE/SE solution(101 
point uniform mesh) of the supersonic nozzle flow- 
field compared with the exact solution represented 
by solid lines. 
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solid line: nonlinear sol. 
dotted line: linear sol. 


(b) outlet 


solid line: nonlinear sol. 
dotted line: linear sol. 


Figure 12: Time history of the acoustic pressure at 
the inlet and outlet for a 101 point uniform mesh. 



dotted line: 101 point mesh 
dash Dot: 201 point mesh 
solid line: exact sol. 




dotted: linear sol. 
dashDot: nonlinear sol. 
A solid line: exact sol. 



Figure 13: The CE/SE nonlinear solutions of the 
acoustic wave pressure p' at the outlet in one time 
period for 101 and 201 point uniform meshes. 


Figure 14: The CE/SE linear solutions of the acous- 
tic wave at t = 40 T using a 101 point uniform mesh 
compared with the corresponding nonlinear solution. 
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